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Matter Molecular Clouds in the Milky Way. 1 Basic 
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ABSTRACT 

We explore some basic observational consequences of assuming that the dark matter 
in the Milky Way consists mainly of molecular clouds, and that cosmic rays can 
penetrate these clouds. In a favoured model of the clouds, this penetration would 
have the following consequences, all of which agree with observation: 

( i) cosmic ray nuclei would be fragmented when they enter a cloud, giving them a 
lifetime in the Galaxy of 10 15 sec (for relativistic nuclei). 

( ii) Pionic 7 rays emitted by the clouds, after pp collisions, would have a flux in the 
galactic plane comparable to the flux from known sources for photon energies >~ 1 
Gev. 

( iii) The heat input into the clouds from cosmic rays would be re-radiated mainly 
in the far infra-red. The resulting radiation background agrees, in both intensity and 
spectrum in different directions, with a known excess in the far i — r background of 
the galaxy over emission by warm dust. 
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1 INTRODUCTION 

Various authors (Pfenniger, Combes & Martinet 1994, De 
Paolis et al 1995, Gerhard & Silk 1996) have proposed 
a viable model in which most of the dark matter in the 
Milky Way is baryonic and consists of a population of self- 
gravitating clouds of mainly molecular hydrogen (possibly 
accompanied by clusters of brown dwarfs). A favoured set 
of parameters for these clouds is: mass M ~ 10~ 3 Mq, ra- 
dius r ~ 10 14 cm, temperature T ~ 10 K, velocity dispersion 
V ~ 200 km • sec -1 , for a halo distribution with a covering 
factor /~5x 10 -4 (Draine 1998). 

One way of testing this model observationally is to 
search for effects arising from the interaction of cosmic rays 
with these clouds. This paper is the first of a series in which 
such an interaction is investigated. In it only the basic impli- 
cations of the interaction are considered. The many details 
involved will be examined in later papers of this series. 

At the outset it has to be decided whether cosmic rays 
can penetrate the clouds. This question is discussed in sec- 
tion 4. If the penetration does occur, a number of observable 
effects would arise for the favoured model of the clouds. In 
this paper the following effects are discussed: 

(i) Cosmic ray nuclei would be completely fragmented 
when they penetrate a cloud. This process may determine 
the survival time of the nuclei in the Galaxy, rather than 
their leaking out of a confinement volume as envisaged in 
the usual picture of cosmic ray propagation. It is straight- 



forward to calculate the mean survival time to, and for the 
favoured cloud model one finds for relativistic nuclei that 
to ~ 10 15 sees. This value agrees well with that derived 
from observations of the radioactive cosmic ray nuclei Be 10 
and Al 26 . By contrast, the survival time in the leakage mod- 
els cannot be calculated a priori. 

(ii) Cosmic ray proton collisions with protons in the clouds 
would lead to the production of neutral pions which would 
rapidly decay into 7 rays. In the favoured model for the 
clouds the resulting integrated 7 ray flux in the galactic 
plane would be comparable to the 7 ray flux due to known 
sources. The clouds would be optically thick for hadronic 
interactions of incident protons with energies exceeding 1 
Gev, and allowance must be made for this in calculating the 
resulting 7 ray flux. A rough estimate of this optical depth 
effect suggests that the excess 7 ray flux in the galactic plane 
would be noticeable only at photon energies exceeding ~ 1 
Gev, and that this excess should be somewhere between a 
factor 1 and 2. This estimate agrees with analyses of the 
EGRET data, which find an unexplained excess in the 7 ray 
flux in the galactic plane of 60% for 7 ray energies exceeding 
1 Gev. 

(iii) Cosmic ray collisions would also heat the clouds, and 
this heat would be mainly re-radiated in the far infra-red 
via molecular line emission. It is convenient to estimate the 
resulting diffuse far i — r radiation F^ r in the galactic plane 
by relating the heat input into a cloud to its emission rate F 7 
for 7 rays with energies exceeding 1 Gev. One obtains in this 



way a conversion factor ~ 0.06 ergs per 7 ray. One can then 
use this conversion factor to derive Fi r from F 7 . We find 
that Fi r ~2x 10~ 5 ergs cm~ 2 sec~ 2 in a sample direction 
in the galactic plane. This intensity, and the spectrum of 
the excess, agree with the properties of a known excess in 
the far i — r emissivity of the Galaxy over that expected 
for a standard warm interstellar dust model. The origin of 
this observed excess is controversial, and we here suggest 
that it is due to the emission from "dark matter" molecular 
clouds powered by cosmic rays. In a later paper in this series 
we study in detail the relation between the intensity and the 
spectrum of the far i — r and 7 ray emissivities of the Galaxy 
as a function of the direction of observation. 

These ideas can also be tested in the future when de- 
tectors, such as Planck, will be able to observe point sources 
with a flux density greater than 1 Jy at 850 pi, even if only 
a few such sources exist over the whole sky. In the present 
model there would be about 10 such cloud sources, with 
correspondingly larger numbers of fainter sources. These 
sources should have S450 ?S Ssso, in contrast to external 
galaxies which have S450 > Ssso- These predictions could 
lead to a decisive test of the model. 

(iv) At high galactic latitudes the averaged column den- 
sity of the clouds is about ten times greater than that of the 
visible interstellar matter. In this case it has already been 
pointed out by De Paolis et al (1999 a, b) that the recently 
discovered 7 ray halo of the Galaxy (Dixon et al 1998) may 
be due to the interaction of cosmic rays with clouds in the 
halo. The implied far i — r emission at high galactic lati- 
tudes then turns out to be of the same order as the observed 
galactic far i — r excess in those directions. 

The plan of the paper is as follows: section 2 is con- 
cerned with some properties of the so-called Extreme Scat- 
tering Events (ESE"s) which are used to determine several 
parameters of the molecular cloud model. In section 3 we 
consider the formation epoch of the clouds which is a par- 
tial guide to their likely metallicity. This metallicity plays a 
crucial role in their far i — r emissivity. Section 4 is devoted 
to a study of ionisation, pionisation and spallation losses of 
the cosmic rays penetrating the clouds, and the value of the 
collision time £0 is determined. In section 5 the temperature 
T of a cloud is derived and also the ratio Fi r /F 7 . The ob- 
served 7 ray emission of the Galaxy and its implications for 
the cloud model are considered in section 6. This is followed 
in section 7 by a discussion of the basic far i — r properties 
of the model which are then compared with observation. 
Finally section 8 contains our conclusions. 

2 EXTREME SCATTERING EVENTS AND 
DARK MATTER MOLECULAR CLOUDS. 

Extreme Scattering Events (ESE"s) (Fiedler et al 1987) con- 
sist of occasional large disturbances to the radio emission of 
point sources. They are believed to be due to plasma ir- 
regularities located in the Galaxy which drift across the ra- 
dio sources. The origin of these irregularities is not known. 
Walker & Wardle (1998) suggested that they are the skins of 
dark matter molecular clouds (hereafter MC's) which have 
been ionised by interstellar ultra-violet radiation. One ad- 
vantage of this suggestion, as Walker & Wardle pointed out, 
is that it leads to a determination of several parameters of 
the MC's, including the crucial result that their total mass 



is comparable to the dark matter mass of the Galaxy (see 
also Draine 1998). 

The ESE's can be used to determine a further parameter 
of the MC's, namely the scale height H of their distribution 
in the Galaxy. This scale height helps to determine the 
number density n of MC's in the galactic plane. This density 
n enters into the calculation of both the collision time to and 
the integrated 7 ray emissivity of the clouds in the galactic 
plane. 

The angular distribution of ESE's in the Galaxy has re- 
cently been considered by Lazio (1999). His statistical anal- 
ysis led him to the conclusion that ESE's are a disk rather 
than a halo population, and so have nothing to do with 
the dark matter problem. However, at my suggestion, he 
kindly considered the intermediate possibility of a flattened 
halo distribution for the ESE's, and in particular the value 3 
kpc for H. His Monte Carlo calculations found (private com- 
munication) that if the radius of the MC distribution in the 
galactic plane is 10 kpc, and H is 3 kpc, then the quadrupole 
moment q of their distribution is within one standard devi- 
ation of the value of q which he derived from the ESE data. 
For a radius of 25 kpc the difference in the g's lies within 
1.5 a. Accordingly we will adopt the value 3 kpc for H. 

This value fits in with a previous analysis of a flattened 
dark halo of our Galaxy (Sciama 1999) where it was shown 
that 3kpc is the smallest allowed value for H (cf. Sack- 
ett 1999). In particular this value is just compatible with 
recent estimates of the total density p of matter near the 
sun, as derived from the vertical motions of nearby stars 
deduced from HIPPARCOS data (Pham 1997, Creze et al 
1998, Holmberg & Flynn 1998). Taking ~ 100 M pc~ 2 for 
the column density of dark matter at the sun on one side of 
the galactic plane (Binney and Tremaine 1987) one would 
have with our adopted value of H a dark matter density 
near the sun of ~ 0.03 Mq pc -3 . Taking I.IM0 pc~ 3 as the 
upper limit of p from the HIPPARCOS data, and 0.04 M 
pc~ 3 for the densities of both stars and visible gas near the 
sun, one just arrives at a consistent set of values. Since each 
MC is assumed to have a mass ~ 1O~ 3 M0, it follows that 
there are ~ 30 MC's pc~ 3 near the sun. 

One may speculate on the origin of the large flattening 
of the MC distribution. One contributing factor could be the 
accretion of interstellar gas subsequent to the formation of 
the MC's. Since the column density of an MC in our adopted 
model is ~ 60 gm cm~ 2 or ~ 10 25 H2 molecules cm~ 2 , which 
is much greater that the column density of the diffuse gas in 
the Galaxy, and since the MC's are moving supersonically 
through this gas, they would simply drill holes in the gas 
and scoop up all the material in their path (Tenorio - Tagle 
et al 1987, Comeron & Torra 1992, Santillan et al 1999). 
Accordingly their column density would increase at the rate 
n g v per unit time, where n g is the number density of the gas. 
Allowing for an increase in n g in the past, and integrating 
over a time of 10 10 years, one arrives at a total accreted 
column density in H2 of ~ 3 x 10 24 cm~ 2 . This is not much 
less than the present column density of an MC, and given 
the various uncertainties in the calculation it is reasonable 
to suppose that the vertical component of the motion of the 
MC's may have been appreciably damped, thereby helping 
to give rise to their flattened distribution. 

Another important consequence of this accretion of in- 
terstellar gas concerns the composition of the MC's. They 



are sometimes considered to have a primordial composition, 
consisting primarily of H and He. However, the presence of 
a significant abundance of heavy elements in the MC's will 
be crucial for our discussion of their far i — r emissivity. As 
we shall see, we will not need a full solar metallicity, but a 
significant contamination could arise from the accreted gas. 
Another source of metallicity will be described in the next 
section. 



3 THE FORMATION EPOCH OF THE MC'S. 

A clue to the likely formation epoch of the MC's can be 
found by noting that if the dark matter in all the galaxies is 
baryonic then the resulting mean density of these baryons at 
the present epoch is close to the value implied by big bang 
nucleosynthesis (BBNS) (Fukugita, Hogan & Peebles 1998). 
Apart from solving the well-known problem "where are the 
baryons?", this result is helpful for determining the forma- 
tion epoch of the MC's because it has recently been realized 
that at a red shift ~ 3 the amount of ionised gas believed 
to be associated with the much smaller amount of neutral 
gas contained in the Lyman a clouds is again close to the 
BBNS value (Rauch & Haehnelt 1995, Giallongo, Fontana 
& Madau 1997). Moreover, the evolution of the Lyman a 
cloud system at lower redshifts is so rapid that the implied 
amount of gas decreases with decreasing red shift on the fast 
timescale of 1 Gyr. Giallongo et al speculate on why the gas 
is disappearing so fast and suggest that it may be driven out 
of the Lyman a clouds by supernova explosions. From our 
point of view it is natural to suggest that one is here ob- 
serving the formation of the MC's. This suggestion would 
imply that most Lyman a clouds are associated with galax- 
ies. This possibility has already been much debated in the 
literature and remains controversial. A recent supporting 
discussion (with references to earlier work) has been given 
by Ortiz-Gil et al (1999). 

One significance for us of these ideas is that at a red 
shift of 3 Lyman a clouds already contained an observable 
abundance of metals (eg. Songaila 1997). This provides a 
second reason for expecting MC's in the Galaxy to contain 
an appreciable abundance of metals. 



4 THE ATTENUATION DEPTH OF MC'S FOR 
COSMIC RAYS. 

If galactic cosmic rays penetrate the MC's, the question 
arises of the attenuation length of the clouds. Let us first 
consider whether this penetration actually occurs. This 
question has already been much discussed, both theoreti- 
cally and observationally, for the known molecular clouds 
in the Galaxy (e.g. Cesarsky & Volk 1978, Lebrun & Paul 
1978, Cesarsky 1980, Berezinskii et al 1990). In their re- 
views of this question Cesarsky (1980) and Berezinskii et al 
(1990) considered that 7 ray observations show conclusively 
that cosmic rays do enter molecular clouds. 

Of course MC's are moving through the interstellar 
medium much faster than the known molecular clouds, and 
it is not clear whether penetration occurs in this case. It 
would be very difficult to carry out a convincing theoretical 
calculation to decide this question a priori and it may be 
better simply to assume that penetration does occur, and 



then compare the implied consequences with observation. 

There are in fact a number of quite different conse- 
quences which can be tested observationally at the present 
time, in particular the lifetime to for cosmic ray nuclei, and 
the 7 ray and far i — r emissivity of the MC's. In addition 
their heating and cooling rates determine their temperature 
T, which is a basic parameter of our adopted model. 

The crucial property of the MC's which controls all 
these phenomena is their column density N. In our adopted 
model N ~ 60gm.cm~ 2 or Nh 2 ~ 10 25 H2 molecules cm~ 2 . 
This is much larger that the grammage required to destroy 
by spallation the cosmic ray nuclei (~ 8 gm.cm~ 2 for light 
and medium mass nuclei and ~ 2 gm.cm -2 for heavy nuclei) . 
It is also comparable to both the stopping power of H2 for 
~ 1 Gev protons due to their ionisation losses, and to the 
attenuation length for the hadronic interactions of Gev pro- 
tons leading to the production of pions (pionisation losses) 
(Review of Particle Physics 1996). Thus on the one hand 
the heat input into the MC's is large. On the other hand, 
self- shielding effects in the MC's will be important. Both 
these facts are used in the present paper, but a full discus- 
sion of them would require extensive calculations which will 
be left for later papers in this series. 

We now consider the calculation of the survival time to , 
and leave for the next section a discussion of the energetics of 
the MC's. As a cosmic ray nucleus propagates in the Galaxy 
it will be scattered by magnetic irregularities in the inter- 
stellar medium. The mean free path associated with this 
scattering process is believed to be about O.lpc (Berezin- 
skii et al 1990). In fact the resulting change of direction 
of a cosmic ray does not matter for the present calculation 
since it is being assumed that the distribution of clouds is 
approximately uniform on the relevant length scale. For a 
relativistic cosmic ray nucleus the mean time to for it to col- 
lide with an MC and be fragmented is given by (7rr 2 nc) _1 . 
For r ~ 10 14 cm and n ~ 30 pc~ 3 one obtains to ~ 10 15 sec. 
This is just the value which has been derived for the mean 
ages of the radioactive cosmic ray nuclei Be 10 and Al 26 with 
energies ~ 1 Gev/nucleon (Simpson & Connell 1998, Web- 
ber & Soutoul 1998). Of course n would be expected to 
decrease with height above the galactic plane, but since the 
storage height of cosmic rays is believed to be less than the 
3kpc scale height of the MC's (Berezinski et al 1990) this 
should be a small effect. We therefore postulate that col- 
lision with MC's provides the main loss process for cosmic 
ray nuclei rather than leakage out of a confinement volume. 

For subrelativistic nuclei of velocity v we would expect 
that to oc 1/v. While the relevant data are not very accurate 
there is some indication that indeed more Be 10 has decayed 
the lower its velocity (Simpson 1983). 

For cosmic ray protons the situation is different. As we 
have seen the column density of an MC is of the same order 
as the attenuation length of H2 for Gev protons. Hence the 
collision of such a proton with an MC would degrade its 
energy but would not destroy it. Allowance should also be 
made for the smaller fractional energy losses of the more 
energetic protons. Again detailed calculations are needed, 
but it is clear that the effective lifetime of cosmic ray protons 
in the Gev energy range could exceed to by a significant 
factor. Moreover much more energetic protons could have 
avoided appreciable energy losses over the whole lifetime of 
the Galaxy. Unfortunately, since the source spectrum of the 



cosmic rays is unknown, it is not possible to test this idea 
observationally. 



5 THE ENERGETIC OF MC'S 

We have already noted that, because MC's in our adopted 
model have a column density ~ 60gm.cm~ 2 , optical depth 
effects would be important for both ionisation and pioni- 
sation losses. Consider first the heating of the MC's by 
cosmic ray interactions. Since 60gm.cm~ 2 of H2 can just 
stop 500 Mev protons due to their ionisation losses (Review 
of Particle Physics 1996) we need to know the energy den- 
sity e CR in cosmic rays of up to about this energy. The 
effective heating rate J per H2 molecule in an MC will 
then be ~ e CR cN~^. For e CR ~ 0.2evcm -3 (Webber 1998) 
and Nh 2 ~ 10 25 cm~ 2 one obtains J ~ 10 _27 ergs.sec _1 H^" 1 . 
There will also be a contribution to J from electrons pro- 
duced by pion and muon decay in the clouds, but for a rough 
estimate this contribution can be neglected. 

To determine the cooling rate one must know the con- 
stitution of the MC's. As discussed previously it will be as- 
sumed that the MC's are contaminated by heavy elements, 
although any dust particles in them must have fallen to their 
centres to avoid opacity effects which have not been observed 
(Gerhard & Silk 1996). The steady state abundances of dif- 
ferent molecules in the MC's represent the outcome of a large 
number of chemical reactions triggered by cosmic ray ioni- 
sation, and have recently been computed for clouds nearly 
as dense as our MC's by Sternberg & Dalgarno (1995) and 
by Neufeld, Lepp & Melnick (1995). The latter authors also 
made extensive calculations of the resulting cooling rate due 
to molecular line emission excited by collisions with H2, al- 
though they used a somewhat small cosmic ray ionisation 
rate. The key point is that, for the large densities and 
optical depths involved in very dense clouds, many of the 
emission lines of a considerable number of molecular species 
become optically thick. An MC is thus rather like a star, 
which is opaque to the radiation produced by its internal 
heat sources. The important consequence for us is that even 
molecular species with a relatively small abundance will con- 
tribute significantly to the emissivity if they are abundant 
enough for their emission to be optically thick. 

Neufeld et al (1995) calculated the emission in all these 
hundreds of optically thick lines, and so derived the cooling 
rate of clouds which are heated and ionised by cosmic rays, 
and are nearly as dense as our MC's. They represented their 
results graphically (for a number of final temperatures) as a 
function of the density nn 2 of H2 molecules in the core of the 
clouds. Gerhard & Silk (1996) represented this cooling rate 



analytically in the form 2.1 x 10 _28 T? 6 /n?i 3 ergs.sec^H" 1 , 
where T = IOT10 K and nu — 10 nn 2 , and used it to de- 
termine the temperature of an MC. For our choice of heat- 
ing rate and MC parameters one obtains T10 ~ 1.9, which 
would satisfy the virial theorem for the clouds (Draine 1998). 
However this calculation should be repeated with a revised 
cooling rate based on the larger cosmic ray ionisation rate 
being used here. 

The strongest emission lines would lie in the far i — r, 
with a cut-off at about 1000 GHz or 300 fi for our choice of 
T. There would be hundreds of lines in this frequency or 
wavelength range, and the resulting spectrum would look 



like a continuum to any detector unable to resolve them. 
To estimate the integrated intensity of this continuum for 
the whole Galaxy one would need to know the spatial dis- 
tribution of both the cosmic rays and the MC's. To avoid 
obtaining a result which is too model dependent we prefer 
to proceed semi-empirically by relating the far i — r emission 
of the MC's to their 7 ray emission. As we shall see, it is 
possible to use observational data on the 7 ray emission of 
the Galaxy as a partial substitute for a spatial model. 

It is pointed out in the next section that the observed 7 
ray emission in the plane of the Galaxy contains an excess 
component for 7 ray energies > 1 Gev. We propose to iden- 
tify this excess as pionic 7 rays emitted by MC's. Now the 
7 ray emissivity of the MC's can be related to their far i — r 
emissivity via a simple conversion factor. To derive this fac- 
tor we note from Mori (1997) that, for a standard cosmic ray 
flux, the number of pionic 7 rays produced with energies > 1 
Gev is 3.2 x 10 _26 sec _1 H 2 " 1 . Because of self-shielding in the 
MC's this flux must be reduced to 1.6 x 10 _26 sec~ 1 Hj 1 (see 
Section 6). Since the far i — r flux ~ 10~ 27 ergs • sec _1 H^" 1 , 
it follows that the conversion factor ~ 0.06 ergs per photon. 
This factor depends on the energy spectrum of the cosmic 
rays but not on their absolute flux. 

6 THE 7 RAY EMISSIVITY OF THE MC'S 

We saw in Section 2 that, in our adopted model, the mean 
density of matter due to the MC's near the sun is nearly 
the same as that of the visible matter in the interstellar 
medium. At first sight this near-equality would lead to an 
approximate doubling of the 7 ray emissivity in the galac- 
tic plane since pionic 7 rays dominate the flux for energies 
exceeding 300 Mev. Yet in their detailed analysis of the 
EGRET data Hunter et al (1987) found good agreement be- 
tween their model, using only visible matter sources, and the 
observational data in the galactic plane, for 7 ray energies 
< 1 Gev. This apparent discrepancy with our model can 
be attributed to self-shielding effects. According to Mori 
(1997) pionic 7 rays of energy < 1 Gev are produced mainly 
by incident protons whose energy lies in the range 1 — 3 Gev. 
These protons are heavily degraded in the MC's, so that a 
substantial fraction of the H2 molecules would be shielded 
from them. On the other hand pionic 7 rays of energy > 1 
Gev are mainly produced by protons with energy > 3 Gev, 
whose fractional energy loss is smaller. Thus shielding ef- 
fects would be less important for these higher energy pro- 
tons. Accordingly some sign of the presence of additional 
higher energy 7 rays would be expected in the observational 
data, if our model is correct. 

In fact, as already mentioned, Hunter et al (1987) found 
a 60% excess in the observed flux in the galactic plane over 
their model calculations for 7 rays with energies > 1 Gev, 
and were unable to give a satisfactory account of this dis- 
crepancy. Other attemps to explain the excess (Mori 1997, 
Strong, Moskalenko & Reiner 1998) have also not been suc- 
cessful. We therefore suggest here that the excess is due to 
pionic 7 rays emitted by MC's powered by cosmic rays. The 
fact that the excess is only 60% rather than closer to a factor 
2 can be taken to indicate some self-shielding even for the 
higher energy protons involved. The corresponding amount 
of self-shielding was used in the previous section to derive 
the conversion factor from 7 rays to far i — r radiation. 



We may add that the mean free path of a 1 Gev 7 ray in 
H2 is 85gm.cm -2 (Review of Particle Physics 1996), so that 
the correction for attenuation of such a 7 ray in an MC is 
not large. However, if we wish to identify the EGRET excess 
with 7 rays from MC's we must impose a strong upper limit 
on the column density r\ of ~ 100 gm cm ~ 2 . On the other 
hand, the absence of an EGRET excess below 1 Gev imposes 
a lower limit on r\ of ~ 30gm.cm~ 2 . Thus rj is somewhat 
constrained. 

According to the recent analysis of Dixon et al (1998) 
the EGRET data also show that our Galaxy posseses a 7 ray 
halo. In this case the relevant column density in the MC's 
would be much greater than that in visible interstellar mat- 
ter, and De Paolis et al (1999 a, b) have already attributed 
this 7 ray flux ( O (10~ 6 — 10~ 5 ) photons cm _2 sec _1 above 
1 Gev) to MC's in the galactic halo. In the next section our 
conversion factor will be applied to this flux as well as to 
the flux in the plane. We shall find that the resulting far 
i — r flux is compatible with observation at high Galactic 
latitudes as well as in the plane. 

7 THE FAR I -R EMISSIVITY OF THE MC'S 

In this section we estimate the far i — r emissivity of the 
MC's and compare the result with observation. To do this 
we use the conversion factor between the far i — r and 7 
ray fluxes from MC's arrived at in the previous section. To 
illustrate the method we estimate the average far i — r flux 
F; r in the range of directions |fej < 2°, £ — 45°. A detailed 
study of various galactic latitudes and longitudes, and of 
the accompanying spectral effects, will be deferred to a later 
paper of this series. 

In our chosen range of directions the average 60% ex- 
cess 7 ray flux ~ 3 x 10~ 4 photons cm _2 sec _ (Hunter 
et al 1987). Accordingly, in these directions Fi r ~ 1.8 x 
10 _5 ergscm^ 2 sec _1 . As already explained, this far i — r ra- 
diation would be concentrated in hundreds of emission lines 
lying mainly in a range dictated by a Rayleigh- Jeans cut- 
off at low frequences and a temperature cut-off at high fre- 
quencies. In our model, with T ~ 20 K, this high frequency 
cut-off would be in the vicinity of 1000 GHz (300 fi). This 
radiation would look like a smooth continuum to a detector 
unable to resolve the closely spaced emission lines. 

Just such a diffuse emission component, cutting off at 
~ 600 GHz, has in fact been discerned in our Galaxy, in the 
form of excess radiation over a warm dust component. This 
dust component has been fitted with a single temperature 
in the range 16 — 21K (for an emissivity law v 2 ). The excess 
component has been detected and interpreted by a num- 
ber of workers (Page, Cheng & Meyer 1990, Meinhold et al 
1993, Reach et al 1995, Lagache et al 1998, 1999, Finkbeiner, 
Davis & Schlegel 1999). For our purposes it is convenient 
to use the analysis of the COBE group (Reach et al 1995), 
whose figure 7 shows clearly the excess component in the 
direction |6| < 3°, t — 45°, as a function of frequency. This 
component has cut-offs at 150 and 600 GHz. Its intensity 
near its peak ~ 1 MJy ster~ , which converts to a total flux 
in this bandwidth ~2x 10~ 5 ergcm~ 2 sec -1 . All these prop- 
erties agree well with those derived above for the far i — r 
emission from MC's in this direction. 

The same test can be carried out at high galactic lat- 
itude, but with less precision. In this case the EGRET 



data have been analysed by Dixon et al (1998) to give for 
the photon flux above 1 Gev the range 10 -6 — 10 -5 pho- 
tons cm _2 sec _1 . As we have seen, this flux has already 
been attributed to MC's by De Paolis et al (1999 a, b). 
The corresponding far i — r flux would be in the range 
6 x 10~ 8 — 6 x 10 _7 ergcm _2 sec _1 . According to Reach et 
al (1995) the excess far i — r component at high latitudes 
~ 0.03 of the value close to the galactic plane, that is, about 
6 x 10 _7 erg cm _2 sec _1 . This value just fits in with the range 
expected from the 7 ray data if MC's are the source of the 
high latitude far i — r excess. 

A further possible observational test of the MC model 
concerns the far i — r emission from a single MC, which to 
many detectors would look like a point source. With 30 
clouds pc~ 3 close to the sun, the nearest MC's would be 
less than a parsec away. In our model an average MC one 
parsec away would have a flux density Ssso at 850 /i ~ 0.5 Jy. 
This would be a bright source for SCUBA, which observes at 
850 /i (Holland et al 1999), but unfortunately this detector 
has a very small field of view. In typical recent work areas 
< lOOarcmin 2 were surveyed (eg. Barger et al 1999 a). In 
such surveys point sources with Sgso in the range 1 — 10 mJy 
were detected (for a review of these observations see Sanders 
1999). Most of these sources are believed to be galaxies with 
red shifts greater than 1, although the identification process 
is a difficult one (Barger et al 1999 b, Richards 1999, Sanders 
1999, Smail et al 1999). The number of these galaxies with 
flux densities Sgso exceeding a given limit is subject to very 
large evolutionary effects, of order 10 3 , as a result of which 
less than 1 percent of the point sources brighter than 1 mJy 
would be MC's. It would be difficult to identify these sources 
with sufficient completeness to pick out this small proportion 
of MC's. However, one way of testing for such a source would 
be to measure its flux density at 450 /i, which SCUBA can 
do. This wavelength is close to the temperature cut-off in the 
MC spectrum, so that for an MC S450 Ss Ssso- By contrast, 
the few sources so far successfully measured at 450 /u,, which 
are presumably galaxies, have S450 ~ 4Ss5o (Blain et al 
1999, Eales et al 1999). 

For the future we note that Planck would easily be able 
to observe the nearby MC's with Ssso ~ 1 Jy, and so test our 
prediction that there should be ~ 10 sources over the whole 
sky with Ssso > 1 Jy- It could also carry out a study of the 
log N — log S relation of these sources and of their fluxes at 
other wavelengths such as 450 /i. 



8 CONCLUSIONS 

If the dark matter in the Galaxy consists mainly of molecular 
clouds (MC's) and if cosmic rays can penetrate these clouds, 
there would be a number of observable consequences. This 
paper has discussed four such consequences for a favoured 
model of the clouds. 

(i) Relativistic cosmic ray nuclei should have a survival 
time in the Galaxy ~ 10 15 sec. 

(ii) There should be an excess in the diffuse 7 ray flux in 
the galactic plane, for 7 ray energies exceeding ~ 1 Gev, of 
order 50% over that due to known sources. 

(iii) There should be an excess in the diffuse far i — r flux 
in the Galaxy over that radiated by warm dust. This excess 
flux should have a cut-off for wavelengths less than about 



450 /i, and should decrease with galactic latitude. Its flux 
density at 850 fi close to the galactic plane at a longitude 
of 45° should be about lMJyster -1 , corresponding to an 
energy flux ~2x 10~ 5 ergcm -2 sec -1 . 

(iv) The MC's would be discrete sources of far i — r radi- 
ation, with about 10 sources in the sky having Ssso > 1 Jy- 
These sources would be expected to exhibit a Euclidean 
log N — log S relation for uniformly distributed sources down 
to a flux level much less than 1 mjy, the present approxi- 
mate limit of the SCUBA counts. At this lower flux level 
less than 1 per cent of the sources at 850 fi should be MC's. 
These MC's would be characterised by having Ssso > S450, 
unlike galaxies which have Ssso ~ 3:8450. 

The first three of these consequences are in agreement 
with observation. The fourth consequence is not contra- 
dicted by the present SCUBA data, and can be tested by 
forthcoming measurements such as those planned for the 
Planck surveyor. 
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